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§1 Introduction 

Let M be an m dimensional compact Riemannian manifold with smooth boundary, 
let V be a smooth vector bundle over M, and let D : C°°(V) -> C°°(V) be an 
operator of Laplace type whose coefficients are independent of the parameter t; 
such an operator is said to be static. There is a canonical connection V on V and 
a canonical cndomorphism E of V so 

(1.1. a) D = -{Tr(V 2 ) + E}. 

Let x — (x\, x m ) be a system of local coordinates on M. We adopt the Einstein 
convention and sum over repeated indices. Fix a local frame for V and expand: 

ds 2 M = g^dx" o dx u and D = -{g^d^ + A^d^ + B) 

where A and B are local sections of TM (3 End(V) and End(V). Let Iy be the 
identity map on V. The connection 1 form u) of V and the endomorphism E 
appearing in equation (1.1. a) are given by 

w 5 = y»s(A v + g^T^Iy) and 

E = B-g"i>{d v u lli + u v w li -u a T vl *)\ 

see [4] for details. Let ';' denote multiple covariant differentiation; we use the Levi- 
Civita connection on M and the connection of equation (1.1. b) determined by D 
to differentiate tensors of all types. If P is a time dependent family of operators 
of Laplace type, then we expand T> in a Taylor series expansion in t to write T> 
invariantly in the form: 

(1.1. C) Vu := Du + J2r>0 t T {Gr, tJ U;ij + ^r^U-i + £ r u] . 
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Typeset by AmS-Te£ 



This setting appears most naturally when denning an adiabatic vacuum in quantum 
field theory in curved spacetime [1]. If the spacetime is slowly varying, then the 
time dependent metric describing the cosmological evolution can be expanded in a 
Taylor series with respect to t. The index r in this situation is then related to the 
adiabatic order. 

Near the boundary, let indices a, b, ... range from 1 through m — 1 and index a local 
orthonormal frame for the boundary; let e m denote the inward unit normal. We 
assume given a decomposition of the boundary dM = Cj^ LJ Cj> as the disjoint 
union of closed sets - we permit CV or Cp to be empty. Let 

(l.l.d) Bu := u\ Cv © (u, m + Su + t(T a u, a + S lU ))\ Cjv 

define the boundary conditions; we can treat both Robin and Dirichlet boundary 
conditions with this formalism. In the following we shall let Bq be the static (i.e. 
time independent) part of the boundary condition; Bqu := u\c?, ® i u ;m + Su)\c^- 
The reason for including a time-dependence in the boundary condition comes e.g. 
from considerations of the dynamical Casimir effect; it takes the form given in 
(l.l.d) for slowly moving boundaries. Here we included only linear powers of t 
because higher orders do not enter into the asymptotic terms we are going to 
calculate. Note that by multiplying B by (1 + T m ) _1 , we can take T m = 0. 
If <f) is the initial temperature distribution, the subsequent temperature distribution 
u,p(t, x) is determined by the equations: 

(l.l.e) (d t + V) U4 ,(t, x) = 0, Bu = 0, and u (O, x) = <j>. 

Let K : <j> — > u^, be the fundamental solution of the heat equation. If V and B arc 
static, then K, = e~ tDe . Let vm be the Riemannian measure on M. There exists a 
smooth endomorphism valued kernel K(t, x, x, V, B) : V s — > V x so 

u<t,(t,x) = (JC4>)(t,x) = j M K{t,x,x,V,B)(t>{x)dvM- 

For fixed t, the operator JC(t) : <f> — > <f>(t, ■) is of trace class. We let 

(l.l.f) a(f,V,B)(t) := Tr L2 (//C(i)) = J M f(x) Tr Vx (K(t, x, x, V, B))dv M - 

The function / 6 C°°(M) is introduced as a localizing or smearing function. As 
t i 0, one can extend the analysis of [6] from the static setting to show that there 
is a complete asymptotic expansion of the form 

(1-1-g) «(/, T>, B){t) ~ E„>o an(/, V, B)t^-™)/\ 

The asymptotic coefficients a n (f,V,B) form the focus of our study. We may de- 
compose a n into an interior and a boundary contribution: 

a n (f,V,B) = a™(f,V)+a d n M (f,V,B). 

The interior invariants vanish if n is odd and do not depend on the boundary 
condition; the boundary invariants are generically non-zero for all n. Let N^(f) 
denote the fi th covariant derivative of the smearing function / with respect to e m . 
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There exist locally computable invari ants a™(x,V) and a™(y,V,B) denned for 
interior points x <E M and boundary points y £ dM so that 

i < = f M f(x) a M(x,V)dv M , and 

a» M (/ ) X>,B) = E„/ SM ^(/)^(l/,©,B)^- 

If 2? and Z3 are static, then these are the heat trace asymptotics which have been 
studied in many contexts previously; a(l,D,B) = Tr^2 e~ tDB . Let Rijki be the 
components of the curvature tensor defined by the Levi-Civita connection and let 
Slij be the components of the curvature endomorphism defined by the auxiliary 
connection V on V. We do not introduce explicit bundle indices for f2jj and E. Let 
L a a be the second fundamental form. Let ':' denote multiple covariant differentia- 
tion with respect to the Levi-Civita connection of the boundary and the connection 
defined by D. We refer to [2] and [4] for the proof of the following result for static 
D; see also related work [3, 7, 8, 9]. 

1.2 Theorem. 

(1) ttf(f,D) = (4tt)-W2 J M fTr(I v ) d u M . 

(2) af{J,D) = (^)- m f 2 \$ M fTY(R l33l Iv + &E)dv M . 

(3) af (/, D) = (47r)- m / 2 ^ J M f Tr{60E. kk + &0R im E + 180E 2 + 30fiyfiy 

+ (12Rijji-kk + SRijjiRkllk — ZRijkiRljkl + ^RijklRijkl)Iv}dvM ■ 

(4) a$ M (f,D,B) = 0. 

(5) af M (f,D,B) = -(4 7 r)( 1 -™)/ 2 i/ CD fTr(I v )dv dM 
+ { ^i- m) ,2ij c j Tr{Iv)dvgM ^ 

(6) a° M (f,D,B) = (4^)- m / 2 i J Cd Tv{2fL aa I v - 3f ;m I v }du dM 
+ (47r)-" l /2 i Tr{f(2L aa I v + 12S) + 2,f, m I v }dv dM - 

(7) af M (/, A3) = -{^) {1 - m)/2 lk fc D Tr{96/£+/(16i?^ 

— 8R a rnma + 7L aa Lbb ~ 10L a bL a b)Iv — 30/ ;m L aa /y + 24/ ;mm 7y }dVQM 

+(47r )(i-m)/2_!_ Tr(96/S + f(i6 RlJJl - 8R amma + UL aa Lbb 
+2L ab L ab )I v + f(96SL aa + 192S 2 ) + f., m (6L aa I v + 96S) 

+ ^f; mm Iv}dVdM- 

(8) af M (f,D,B) = (4^)- m /2_i_ Tr{f(-120E. m + l20EL aa ) 

f ( ^-&Rijji;m ~t~ ^^Rijji-^aa ~t~ ^-Ramam-^bb ^^Rarabra^ab ~t~ ^Rabcb-^ac 
-\-2AL aa -bb + 2iL aa LbbL cc Y^abLabLcc + -55- L a bLb c L ac )Iv — 18>0f; m E 

+.f;m(—30Rijji — ^Y-L aa L b b + ^-L ab L ab )I v + 2Af. mm L aa I v 

— 30f;iimIv}dl'dM 

+ (4^)-™/2_i_ Tr{/(240£ ;m + 120EL aa ) + f(42R ijji;m + 24L aa:bb 

-\-10Rij jiL aa -j~ ^RamamLbb ~~ 12i? am fr m Z/ a 5 H~ ^Rabcb-^ac H 3" ^aa-^bb-^cc 

-\~SL a bL a bL 

cc ~t~ g LabLbc-^ac )I V + f(720SE + WOSRijji + U4SL aa L bb 
+A8SL ab L ab + mS 2 L aa + 480S 3 + 120S :aa ) + / ;ro (180£ + 72SL aa 
+240S 2 ) + /:,,,; 30/.', ; „ + \2L aa L bb + l2L ab L ab )I v + 120 f. mm S 

+ 24:f;mmLaaIv + 30/ ; jj m Jy }dv dM . 



The main result of this paper is the following result which extends Theorem 1.2 to 
the time dependent setting: 
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1.3 Theorem. 



(1) a™(f,V)= a M(f,D). 

(2) a™(f,V) =a 2 M (/,i)) + (4 7 r)-™/ 2 i/ M /Tr(|g Mi )^M. 

(3) af{f,V) = af{f,D) + (4,)-/^/ M / Tr(f g 1M g ltjj + f Mj Mj 
+60^ 2 ,ii - 180£i + l5Qi^iRjkk 3 - SOGi^jRikkj + 90Gi,uE + 60^1^^ 
155i, u-jj — 30Gi.ij;ij)dvM ■ 

(4) a d n M {f,V,B) = a d n M (f,D,B ) for n < 2. 

(5) a d 3 M (f,V,B) = af M (f,D, fi ) + (4tt)( 1 -'™)/ 2 ^ J Cd f Tr(-240 liaa )^ aM 

+ ^ n )^-rn)/2^_J c j Tr{24 g iaa)dudM _ 

(6) 4 M (f, V, B) = a d 4 M (f, D, B ) + (47t)-™/2_!_ J Cd Tr{/(3O0 liOa L 66 
-6001 .abL a b + 3O0i 305i ,aa;m 

+ O0i 

,am;a 

+ /;m(-455l,aa + ^Gl,mm)}dv aM 

+ (^y m/2 m> fc N Tr{/(30ei,aai b6 + 120g hmm L bb - 15O0i, a6 L a6 

,mm;m 

+ 6O0i 

,aa;m 

+ O0i 

.am; a 

+ 150JPi im + 18OS01 

.aa 

-18O50i, mro + 360Si + 0T a:a ) + / ;m (450i iOO - 45Gi, mm )}dv dM ■ 

Here is a brief outline to this paper. In §2, we use invariance theory and dimensional 
analysis to study the general form of the invariants a n (f,V,B). We shall use B~ 
for Dirichlet and B + for Robin boundary conditions. We shall show, for example, 
that there exist constants c and ef so that: 

af (/, V) = af (/, D) + (47r)-™/ 2 i J M f Tr(c Gi,u)dv M and 
a! M (f,V,B) = a! M (f,D,B ) + (4tt)-(™- 1 )/ 2 ^ / Cd / Tr(e^Gi, aa )du dM 

we refer to Lemma 2.1 for further details. The interior invariants will be described 
by constants {ci}^£ , the boundary invariants for Neumann boundary conditions 
will be described by constants {ef}}^, and the boundary invariants for Dirichlet 
boundary conditions will be described by constants {e7}jii- We use the localizing 
function / to decouple the interior and the boundary integrals; with the exception 
of Lemma 2.4, there is no interaction between the unknown constants {ci}, {e7}, 
and {e + }. A priori, those constants could depend on the dimension. In Lemma 
2.3, we will use product formulas to dimension shift and show the constants are 
dimension free. We complete the proof of Theorem 1 .3 by evaluating these unknown 
constants; the values we shall derive are summarized in Table 2.2. 
We use various functorial properties to derive relations among these constants. For 
example, in Lemma 2.4, we use the product formulas of Lemma 2.3 to show that 
C5 = 10cq. The functorial properties that these time dependent invariants satisfy 
and which are discussed in §3-§6 are new and have not been used previously in 
other calculations of the heat trace asymptotics. Thus we believe they are of in- 
terest in their own right. It is one of the features of the functorial method that 
one has to work in great generality even if one is only interested in special cases. 
We found it necessary, for example, to consider the very general time dependent 
boundary conditions of equation (1.1. d) to ensure that the class of boundary con- 
ditions was invariant under the gauge and coordinate transformations employed in 
§4 and §5. We work with scalar operators as the (possible) non-commutativity of 
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the endomorphisms in the vector valued case plays no role in the evaluation of a n 
for n < 4. 

We summarize the five functorial properties we shall use as follows. In §2, we 
consider a product manifold M — Mi x M 2 where dM 2 is empty, and an operator 
of the form T> = T>\ <g> 1 + 1 <g> T> 2 . In Lemma 2.4, we show that 

a n {hf 2 ,V,B) = J2 p+q=n a p (f 1 ,V 1 ,B)a q (f 2 ,V 2 ). 

In §3, we rescale the time parameter t. Let D and B be static operators. Let 
V := (1 + 2at + 3/3t 2 )D. In Lemma 3.1, we show that: 

a 2 (f, V, B) = a 2 (f, D, B) - faa (f, D, B) 
a 3 (f, V, B) = a 3 (f, D, B) - ^aa, (/, D, B) 

a 4 (f, V, B) = a 4 (f, D, B) - ^aa 2 (f, D, B) + C-^^a 2 - f (3)a (f, D, B). 

In §4, we make a time dependent gauge transformation. We assume D and B 
are static. Let V e :— e~* e *De te * + g^. We also gauge transform the boundary 
condition B to define B e . In Lemma 4.1, we show that : 

d-{a n (f,V e ,B e )}\ e=0 = -a n _2(/*,£>,B). 

In §5, we make a time dependent coordinate transformation. Let A be the scalar 
Laplacian and let B be static. Let $ e : {t,x\,x 2 ) — > (i, x\ + tgR,x 2 ) where g is 
an auxiliary parameter. We set V B := $*(<9 t + A) - d t and B e := &* g (B). Let 
dvM := gdx 1 dx 2 . In Lemma 5.1, we show that: 

f e {a n {f,V e ,B e )}\ e = - -ia n _ 2 Gr 1 0iG//H),A,B). 

In §6, we assume given a second order operator Q which commutes with a static 
operator D of Laplace type. We define D g := D+gQ and define a suitable boundary 
condition B e . We also define V e := D + 2tgQ and show 

f e {a n {f,V ei B)}\ e=Q = f e {a n _ 2 (f,D e ,B e )}\ e=0 . 

In each section, we use the relevant functorial properties to derive relations among 
the unknown coefficients; these relations are contained in Lemmas 2.4, 3.2, 4.2, 
and 5.2. These relations suffice to determine the unknown coefficients and thereby 
complete the proof of Theorem 1.3. As the computations are somewhat long and 
technical, we have derived more equations than are needed as a consistency check; 
this is typical in such computations. 

§2 Invariance Theory, dimensional analysis, and dimension shifting 

We begin the proof of Theorem 1.3 by establishing the general form of the invariants 
and a® M for n < 4. Let (D,Bo) be the static operator and boundary condition 
determined by (V,B). 
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2.1 Lemma. There exist constants so that 



(1) a M (.f, V) = a M (f, D) and af M (f, V, B) = af M (f, D, B ) for i < 2. 

(2) a™(f,V) =aM{f,D) + {l-K)-™/*\j M fTr{c Q Q 1 , ii }dv M . 

(3) of (/, V) = af{f, D) + (4ir)- m / 2 ^ f M f Tr {aG^uGuj + <^J..,Ah ,,, 

+ C3G2M + C4(?i + C$Gl,iiRjkkj + CeGl,ijRikkj + CjGl.iiE + C%T\^i 

+ c 9Gi,u-jj + c w Gi,ij;ij}dvM ■ 

(4) a d 3 M (f,V,B) = al M (f,D,B ) + (^) (1 - m)/2 ^ J Cd fTr(e^Gi,aa 
+e^Gi, mm )dy aM + (47r) (1_m)/2 3g4 J Cn /Tr(ef C/i >00 + e\G\, mm )dv dM ■ 

(5) 4 M {f,V,B) = a 4 (f,D,B ) + (47r)- m /2_i_ J Cd Tv{f{e^G haa L bb 

+ e 4 Gl,mmLbb + 65 Gl,abL a b + e g Gl,mm;m + e^ Gl.aa;m 

,am;a ~\~ 6g •T'l.m) "I" f;m{^l()Gl,aa + enGl,m.m)}dVdM 

+ (4^)- m /2 _i_ Tr{/(e+ai,aaib6 + eJei, mm L 66 + e+£i, o() L a6 

,mm ~t~ e i4*5'l ~t~ e i5^a:a) H~ /;m( 

Proof. We use dimensional analysis - this involves studying the behavior of these 
invariants under rescaling and is described in [4] in the static setting. We assign 
weight 2 to _R, f2, S and T a and weight 3 to Si. We assign weight 1 to S and L ab . 
We increase the weight by 1 for each explicit covariant derivative which appears. 
Thus, for example, the terms E-kk, SlijSlij, and RijkiRijki are all of degree 4. The 
integrands appearing in and a^ M are weighted homogeneous of degree n and 
n— 1. The structure groups are O(m) and 0(m— 1) respectively. H. Weyl's Theorem 
[10] shows that all orthogonal invariants are given by contractions of indices. The 
assertions of the Lemma now follow by writing down a spanning set for the space 
of invariants. We remark that since Gi,ij — Giji, the invariant Gi,ij^ij does not 
appear. □ 

We will complete the proof of Theorem 1.3 by evaluating the unknown coefficients 
of Lemma 2.1. The remainder of this paper is devoted to deriving the values in the 
following table: 



Table 2.2 



co = | 


Cl = ^ 

t-l 4 


C2 = f 


c 3 = 60 


c 4 = -180 


c 5 = 15 


c 6 = -30 


c 7 = 90 


c 8 = 60 


Cg = 15 


cio = -30 




= -24 


e 2 - = 


eg = 30 


el = -60 


e^ = 30 


e 6 " = 30 


ef = -30 


eg =0 


eg = -30 


efo = -45 


en = 45 




e+ = 24 


e+ = 


e+ = 30 


e+ = 120 


e+ = -150 


e+ = -60 


e^~ = 60 


eg = 


ej = 150 


e+ = 45 


e+ = -45 


ef 2 = 180 


e+ = -180 


e+ = 360 


et,=0 









The (possible) non-commutativity of the endomorphisms in the vector valued case 
plays no role in the invariants of Lemma 2.1. We therefore suppose V to be the 
trivial bundle hence forth and omit the trace from our formulas to simplify the 
notation as we will be dealing with scalar operators on C°°(M). We also set e^ = 
for i > 12 to have a common formalism; these constants describe invariants which 
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involve S, Si, and T a and which are therefore not relevant for Dirichlct boundary 
conditions. 

A-priori, the constants Cj and ef might depend upon the dimension. Fortunately, 
this turns out not to be the case; the dependence upon the dimension is contained in 
the multiplicative normalizing factors of (47r)*. Let T>i be smooth time dependent 
families of operators of Laplace type over manifolds Mj for i = 1, 2. We suppose 
M 2 is closed. Let M := M\ x M 2 , let T> := T>i+T> 2 , and let the boundary condition 
for M be induced from the corresponding boundary condition for Mi. 

2.3 Lemma. Adopt the notation established above. 

(1) a™{hf 2 ,V) = E P+q=n < 1 (.h^iKH.f2,v 2 ) 

(2) a d n M {hh,V,B) = E P+q =n4 Ml (-h^ B K 2 (fM- 

(3) The constants of Lemma 2.1 do not depend upon the dimension m. 

Proof. We use equation (1.1. e) to check that u < p 1 . ( f >2 — u ( p 1 ■ u^. This shows the 
kernel function on M is the product of the corresponding kernel functions on Mi 
and on M 2 ; assertions (1) and (2) now follow. Let (M ,T>m ,B) be given. Let S 1 
be the unit circle with the usual flat metric and usual periodic parameter 9. Let 
D s = —dg on the trivial line bundle. Let T> MxS i = V M + D S . Then a p (9, D s ) = 
for p > and a a (6,Ds) — (47r) -1 / 2 ; see [4] for details. Thus p — n and q = in 
assertions (1) and (2) so a n (fi,V MxS i) = (47r) _1 / 2 a„(/i, £>m, It now follows 
that Ci(m + 1) = Ci(m) and ef(m + 1) = ef(m). □ 

We use the product formulas of Lemma 2.3 to prove the following Lemma: 

2.4 Lemma. We have c\ = 5cq ; c 5 = 10c , c-j — 60co, = — 16c , ef = 20co, 
ef = — 30c , e\ = 16co, ef = 20co, ef = 30c 0; and ef 2 = 120c . 

Proof. We apply Lemma 2.3 and study the cross terms arising in a p + q (/i/2, T>, B) 
from a p (fi, T>i, Bi)a q (f2, T> 2 )- We let indices r and s index Mi and indices u and 
v index M 2 . We use Theorem 1.2 and equate coefficients of suitable expressions to 
derive the following systems of equations from which the Lemma will follow: 



2ci - 360(ic )(ic ) [ftf2g1.rrG1.uu] 


c 5 = 360(i)(ic ) [J1f2Rrs.srQ1.uu] 


c 7 = 360(ic ) [fifiEiQi^u] 


e±= 384(±i)(i Co ) [/i/ 2 ^i,„„] 


eg = 360(i)(|c ) [f1f2L.rrQ1.uu] 


ef = 360(±|)(ic ) [/i ;m / 2 ^i, u „] 


e+ = 360(2)(ic ) [fSGi.uu] □ 





§3 Rescaling the time parameter 

Let D and B be static. Let a, (i G M. We define a time dependent family of 
operators of Laplace type by setting: T> := (1 + 2crf + 3[3t 2 )D. 

3.1 Lemma. 

(1) a 2 (f, V, B) = a 2 (f, D, B) - faa (f, D, B). 

(2) a 3 (f, V, B) - a 3 (f, D, B) - ^aai(f, D, B). 

(3) a 4 (/, V, B) = a 4 (/, D, B) - ^aa 2 (f D, B) + { r -^±^o? - f /3)a (/, D, B). 
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Proof. Let u = e- tDB cj) and let u(t, x) := u (t + at 2 + /3t 3 , x). Then: 

Vu(t, x) = (1 + 2at + 3f3t 2 )(Du a )(t + at 2 + /3t 3 , x) 
d t u(t,x) = {l + 2at + 3l3t 2 ){d t u Q )(t + at 2 + f3t 3 ,x). 

This shows that (d t + V)u = 0. Since ^(0,2;) = u (0, x) — <p(x) and Bu = 0, the 
relations of equation (1.1. e) are satisfied so 

K (t, x, x, V, B) = K(t + at 2 + (3t 3 , x, x, D, B). 

The Lemma will then follow from the expansions: 

a(f, V, B){t) ~ £„ t- m l 2 {\ + at + pt 2 )^- m )l 2 a n {f \ D, B)t n / 2 
(1 + at + /3t 2 y ~ 1 + ajt + + j/3)t 2 + 0{t 3 ) □ 

We apply Theorem 1.2 and Lemma 3.1 to derive the following relationships: 
3.2 Lemma. 

(1) c = |, ci = f , c 2 = f , c 3 = 60, c 4 = -180, c 5 = 15, c 6 = -30, c 7 = 90. 

(2) ef = ±24, ef = 0, ef = 30, ef + e± = -30, e± = ±45, e± = ^45. 

(3) e+ = 180, e+ - -180. 

Proof. We have C/i^j = —2agij, T\.i = 0, ^2,ij = —3f3gij, and £i = — 2aE. Thus 
Gi,u;jj = 0, Gi,ij-,ij — 0, and .Fi,^ = 0. We equate coefficients of suitable ex- 
pressions in Lemma 3.1 to derive the following systems of equations from which 
the Lemma will follow. Note that since m is arbitrary, equations involving this 
parameter can give rise to more than one relation. 



-2mc = -6f 


[af] in af 


4(m 2 Cl + mc 2 ) = 360 m( " g l+2) 


[a 2 /] in af 


-3mc 3 = -360f 


[/?/] in af 


-2(c 4 + mc 7 ) = -360^6 


in af 


-2(mc 5 + c 6 ) = -360%^ 


[a/fiijji] in af 


-2{(m - l)ef + e±} = -384(^)(±|) 


[a/] in af M 


-2{(m - l)e± + e± + ef} = -360(^)(|) 


[afLaa] in af M 


-2{(m - l)ef + 4} = -360(^)(±I) 


[a/ ;m ] in af M 


-2{(m-l)e++e+} = -360(2^)(2). 


[afS] in af M □ 



§4 Time dependent gauge transformations 

Let V e := e - t<? *L>e t<? * + a*. If S u — ii;m ± 'S'a is the Robin boundary operator, 
we gauge transform the boundary condition to define B e := V m + S + tSi with 
Si = Q^-m'i the Dirichlet boundary operator is unchanged. 
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4.1 Lemma. We have -^{a n (f,V e ,B e )}\ e=0 = -a„_ 2 (/* ,D,B). 

Proof. Let u := e~ tDe (j) and let u :— e~ te *u - We show u satisfies the relations of 
(1.1. e) by computing: 

d t u(t,x) = e- te * {d t - g^)u , V e u(t,x) = e- te *(D + Q*)u , 

{dt + V e )u = e- te *{d t + D)u = 0, and u(0, x) = u (x) = 4>{x). 

Dirichlct boundary conditions are preserved. With Robin boundary conditions, 
u. m + Su + tSiu = e~ te *(u 0; m - tg^. m u + Su + tg^. m u ) = 0. 

Thus K(-,T> g , B e ) = e-^K^D, B). The Lemma now follows. □ 

We use Lemma 4.1 to obtain some additional relationships: 

4.2 Lemma. We have c 8 = 60, = —30, and — 2e$ = 60. 

Proof. Let * vanish on dM. We apply Lemma 4.1 with M = [0, 1] and D = -d%. 
We work modulo terms which are 0(g 2 ) and compute: 

D e = D + g^- 2tg^, d g - tg^.gg, 

B+ = V m + S + tg<f. m , 5i = e* ifl , 

E = T x . m = -2g^ :8 , £ x = -g^-ee- 



We study §-{af }| e=0 and f {a d 4 M }\ e=0 : 



f e {60E, u }\ e=0 = -60<i> 


^{-180fi}| e=0 = ISO^.gg 






^{(-120-, 240+)£ ;m }| e=0 = (120- , -240+)* ;6 





Here the notation (—120 , 240 + ) indicates that the coefficient for Dirichlct B and 
Neumann B + boundary conditions is —120 and 240. As — a^(f^ } D) = and 

-af^/*, = -3M( 47r ) _1/2 /a M ±180 (/*);™' we use Lemma 4.1 to derive 

the following equations from which the Lemma will follow: 

= -60+ 180 - 2c 8 , 

-180= -2e+ + e+ - 240, 

180 = 120 - 2eg. □ 

§5 Time dependent coordinate transformations 

In this section, we study time dependent coordinate transformations and make a 
coordinate transformation that mixes up the spatial and the temporal coordinates. 
This technique was also used in [5] to study the heat content asymptotics. We work 
in a very specific context but note the Lemma holds true in much greater generality. 
Let M := S 1 x [0,1] with ds 2 = e 2 '^dx\ + e 2 ^dx\. Let dv M := gdx 1 dx 2 . Let 
5 € C°°(M) have compact support near some point P £ M. Let A be the scalar 
Laplacian and let B be a static boundary condition. Define: 

$ e (t, xi,x 2 ) := {t, xi + tgE, x 2 ), 

V e := $* e (dt + A) - d t , and B e := $*(£). 
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5.1 Lemma. We have -§^a n (f,T> e ,B e )\ e = = -\a n -2{g 1 d 1 (gfE),A,B). 

Proof. Let u(t,xi,x 2 ) := {&*(e~ tAB <f))}(xi, x 2 ). By naturality, u satisfies the rela- 
tions of (1.1. e). As the static operator determined by T> e is A+ lower order terms, 
dvM is independent of g. Thus 

K(t,xi,x 2 ,x 1 ,x 2 ,V e ,B e ) = K(t,xi + gtE(x 1 ,x 2 ),x 2 ,x 1 ,x 2 ,A,B). 

We set X\ — x\ and x 2 — x 2 . We work modulo terms which are 0(g 2 ) and expand 
in a Taylor series to compute: 

a(f, T> g , B e )(t) = J M f(xi,x 2 )K(t, xi,x 2 , xi,x 2 ,V g , B e )dv M 

= Im f{xi,x 2 )K(t,xi + gtE,x 2 ,xi,x 2 ,A,B)gdxidx 2 

= Jm^^ 1 ' x 2 )K(t,xi, x 2 , xi,x 2 , A, B) 

+ tgfZdiK(t,x ll x 2l yi,x 2l A 1 B)\ Xl=yi }gdxidx 2 . 

As Ab is self adjoint, the heat kernel is symmetric. Thus we have: 

a(f,V e ,B e )(t) = f M {f(x 1 ,x 2 )K(t,x 1 ,x 2 ,x 1 ,x 2 ,A,B) 

+ \tgfEd\K{t, xi,x 2 ,Xi,x 2 ,A, B)}gdx\dx 2 

= I M {f(xi,x 2 )K(t,xi,x 2 ,xi,x 2 ,A,B) 

- ^tgg~ 1 d 1 (gfE)K(t, xi,x 2 ,xi,x 2 , A, B)}dv M 

= a(f,A,B)(t) - Itgaig-^igfE), A,B)(t). □ 

We use Lemma 5.1 to complete the proof of Theorem 1.3 by completing the calcu- 
lation of the coefficients Ci and ef . 

5.2 Lemma. 

(1) eg = 15 and c w = —30. 

(2) = —60, e r = 30, = 30, e r = —30, and e$ = 0. 

(3) e\ = 120, e+ = -150, e+ = -60, e\ = 60, e+ = 0, e+ = 150, e+ = 360, 
and ef 5 = 0. 

Proof. We introduce an auxiliary parameter e and work modulo terms which are 
0(e 2 ) + 0(g 2 ). Let 

ds 2 :=e 2 ^dx 2 + e 2 ^dx 2 . 

The Laplacian A = —g^digg^dj can then be expressed in the form 

A = -{e- 2 ^0 2 + e- 2 ^d 2 + e(i> 2/1 - + e{4> 1/2 - </> 2 / 2 )5 2 }. 

Let x\, x 2 ) = (t,xi + gtE,x 2 ). Let En = diE etc. As $ e is a diffeomorphism, 
we can pull back both differential forms and differential operators. We compute: 

= d t - tgZ/rfu $*(9 2 ) ee 9 2 - isS/aft,**^) ee 9 t - gEd^ 

The operator £> e := <&*(<9 t + A) - <9 t is given by: 

V e =A + tg{e- 2 ^ [2E fl d 2 + + e^^E^d^ + S /22 di]} 

+ tge{2^p 1/1 Ed 2 + 2ip 2/1 Ed 2 + E /1 (ip 2/1 - ^i/i)^ 

- S(^2/n - ^i/n)di + E /2 (ip 1/2 - ip 2/2 )di 

- S(^i/i 2 - ^2/12)^2}. 

The tensors £?, and £1 are therefore given by: 
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D = A - Q~di 


uj? = \e^Q~ 

1 2 [ - 


gi u = e -Wi2 g Z /1 + 2eip 1/1 Q~ 




Q x 22 = 2eil> 2/ie E 


= e- 2 ^ Q Z /2 


E = -\qB./i - 2 -e(ip 1/1 + ip2/i)e^ 


£ 1 =0 



To compute J 7 , we must express partial differentiation in terms of covariant differ- 
entiation. Since ui is linear in g, it plays no role. The Christoffel symbols of the 
metric, however, play a crucial role. We compute: 

Qi n /;ii = (Gi U df - 2gE /1 e^ 1/1 d 1 + 2gE /1 e^ 1/2 d 2 )f 
20i, 12 / ; i2 = (25i, 12 a 1 a 2 - 2gZ /2 ei> 1/2 d 1 - 2gZ /2 e^ 2/1 d 2 )f 
Qi 22 f,22 = Gi 22 d 2 2 f 

We use this computation to determine the tensor T\: 

+££{002/1 - ^i/i)2/i - (-02/11 - -0i/n)2 

+ (V>l/2 - "02/2)2/2 + 2^1/lS/i + 2^1/ 2 S /2 } 
T X 2 = £^{-(^1/12 - "02/12)2 - 2^1/2^/1 + 2^2/lS/ 2 } 

We now prove assertion (1). Let P e int(M). Let e4>i(P ) = eip 2 (P) = 0. We 
study monomials S/ m and ^2/111 ^ appearing in ^{^(Oll^o- Let 1Z — E or let 
1Z = Rijji- We integrate by parts to define A[R] by the identity: 

_ T2 Im 9~ 1 d 1 {gf'S)'Rdu M = J M fA[1l}dv M ; then 

- XGr^Gz/S), A) = (47T)- 1 ^ J M /,l[625 + P im \dv M . 

We have -R^ji = — 2e^ 2 /n + We compute: 



^{60E. u }\ e=0 = 


-30S /m 


-3O£-02/lllS 


+ 


■^{Q0^Fl,i;i}\g=0 = 


60S /m 


-60£V2/lll2 


+ 


■^{c9Gl,ii;jj} e=0 = 


2cgS/ m 


+2c 9 e-02/iiiS 


+ 


^{ciO01,y;y}| e =O = 


2cioS/m 


+Oci O £-0 2 / m S 


+ 


.a[6.e] = 


0H /m 


+Oe-0 2 /iiiS 


+ 




OH/m 


-60£V2/ni5 


+ 



We use Lemma 5.f to relate the coefficients of /S/ m and ,f02/iii£ and establish 
the following relationships from which assertion (1) follows: 

-30 + 60 + 2c 9 + 2ci = and - 30 - 60 + 2c 9 = -60. 

We now study the boundary terms. We pullback the Robin boundary operator 

$* e (e- e ^d 2 + S) = e~ £ ^ teB {B - e-^HgZ^d! + tg~(Se^ 2/1 + S fl )} 
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to determine the tensors 

T 1 = -e-^ 2 Q~ /2 and S 1 = gZ{e^ 2/1 S + S fl ). 

We have in = —eipi/ 2 - We study the terms comprising J^{af M (f,T> g , B g )}\ e= o- 
At the point of the boundary in question, we suppose etpi(P) = eip2(P) = 0. 
^{(-120-,240+)/S ;ro }| e=0 

= (60", -120+)/{H /12 + H1/12 + «/>2/i 2 )S + (e^i/i + ^2/1)^/2}, 
£{120 fEL aa }\ e=0 = 60£M /2 S /1; 
^{720/5^} I e=0 = -360/5{Sj + e(^ vi + ^/i)S}, 
^{e^/^i,aaife6}| e =o = e^/(2S /1 )(-£^i/ 2 ), 

J^{ e 4 fQl,rnmLbb\\g=0 = 0, 

^{4/Si,-6i-6}l«?=o = e±/(2S /1 )(-£Vi/ 2 ), 

^{e^/^i, mm;m }| e =o = eff(2etp 2 /i2^ + 4£^ 2 /iS/ 2 ), 

^{ e 7 /^i,oo;m}| e =o = eff{2E /12 + 2£^i/i 2 S + 2eV>i/iS /2 - 2£-0 2 /iS /2 }, 

^{eg /Si,om;o}| e =o = eJ/{-£-0 2 /i2/2 + S /12 + eVi/i s /2 ~ 2e^ 1/2 S/i}, 

^{ e 9 /• ?7 l,»n}le=0 = eg /{-(e^l/12 - £-02/12)2 - 2^1/25/! + 2£-0 2 /lS/ 2 }, 

^{e+ fSg haa }\ e=0 = e+ 2 f{2E n S + 2e^ 1/1 SS}, 
f e {et 3 fSg hmm }\ e=0 = e+ 3 f{2e^ 2/1 ES}, 
£{etJSi}\ e =o = e+/S{£V 2 /iS + S fl }, 
£{et 5 f T a:a}\ e =o = ef 6 f (etp 2/ - H /12 - etp in E /2 ), 
^{(±180)/ ;mJ B}| e=0 = T90/ ;m {S A + (e^ in + £^ 2 /i)S}, 
^{ef /;m^i,aa}|e=o = e%f. m (2E /1 + 2e^i/iS), 
^{e n / ;m ^i, mm }| e =o = e 11 / ;m 2£'02/i2. 

We must also study the boundary terms comprising — \a 2 M {-)- As when studying 
a^ 1 , we integrate by parts to define *4 and compute: 

A[2fL aa ] = -60eM/i 2 S, 

^[12/5] = -360{S £ /^2/i ~ /SS/i}, 

-4[±3/ ;m ] = T90{(£^ 1/12 + £^ 2 /i 2 )/S + 2£0 2/1 (/ ;m S + /S /2 )}. 
We established the following relations in Lemmas 3.2 and 4.2: 

= 30, + ef = -30, e+ - 2e+ = 60 and = -30. 
We use Lemma 5.1 to derive the following equations and complete the proof: 



(60-,-120+)+4e^ 


" 2e 7 ± - 


- + 2e* + 


ef B = T180 


[/£-0 2/ lS /2 ] 


(60-,-120+) + 2e^- 




T90 




[M2/12=] 


(60", -120+) + 2ef - 


± 

" e 9 = 


-60^90 




[/^i/i 2 S] 


(60", -120+) + 2ef - 


Fef - 


ef 5 = 




[/S/12] 


2c^ 2cg — 





[/#l/2S/l] 


e+ = 360 




-360 + 2e£j + e+ = 


-360 


[/e^/iSSJ 


-360 + 2e^ = 


= [/s A s] 


T90 + 2efi = T180 




[/;m£02/l2] 


T90 + 2ef = 


[/jmS/x] 
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§6 Commuting operators 

We conclude this paper by deriving a final functorial property. The equations which 
can be derived using this property are compatible with the values for the constants 
d and ef previously computed; they are omitted in the interests of brevity. 

6.1 Lemma. Let D be a self-adjoint static operator of Laplace type and let B 
be a static boundary condition. Let Q be an auxiliary self-adjoint static partial 
differential operator of order at most 2 which commutes with D and with B. Then: 

§- e {a n {f, D + 2tgQ, B)}\ e=0 = £{a„_ 2 (/, D + gQ, B)}\ e=0 . 

Remark. If we take D = Q, then D(g) = (1 + 2tg)D. By Lemma 3.1, 

f g {a 4 (f, (1 + 2tg)D, B)}\ e = = ^a 2 (f, D, B). 

On the other hand, clearly a n (f, (1 + g)D,B) = (1 + g) (n " m)/2 a„(/, D, B). Thus 
we may show that Lemma 6.1 is compatible with Lemma 3.1 in this special case by 
computing: 

(1 + Q)D, B)}\ e=0 = ^a 2 (f, D, B) = ^{a 4 (/, D + 2tgD, B)}| e=0 . 

Proof. Let Ki{i) := (1 - t 2 gQ) e - tDB . Then /Ci(0) is the identity operator and: 

{d t +D + 2tgQ)(l-t 2 gQ)e- tDe 
={-2tgQ - (1 - t 2 gQ)D + D(l - t 2 gQ) + 2tgQ{\ - t 2 gQ)} e - tDB 
= -2t 3 g 2 Q 2 e- tD *. 

There exists a constant C and an integer \x so that we have the estimate in a suitable 
operator norm: 

\-2t 3 g 2 Q 2 e- tDs \<Ct-"g 2 . 

Thus since we are interested in the linear terms in g, we may replace the fun- 
damental solution of the heat equation K.(t) for D + 2tgQ by the approximation 
(1 — gt 2 Q)e~ tL>B . There is an asymptotic expansion of the form [4]: 

Tr L2 (/Qe- tDB ) ~ £„>o t("- ro " 2 )/ 2 a„(/, Q, D,B). 

We equate coefficients of ^f™-" 1 )/ 2 i n the asymptotic expansions to see 

£ e {a n (f,D + 2tgQ,B)}\ e=0 = -a n - 2 (f,Q,D,B). 

Since Q and D commute and since Q and B commute, we complete the proof by 
computing: 

Zn> £ e Mf,D + gQ,B)}\ e = t("-™V 2 ~ 9_ {T r L2 (fe-^ D+ ^M)}\ e =o 

= Tr L 2(-tfQe- tD *) ~ - E„>o a ™(/> <5. D , S)^"" m )/ 2 so 
£{a n (f, D + gQ, B)}\ e=0 - -a n (f, Q, D, B). □ 
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